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Abstract 

A method for Including elastic effects In 
steady, transonic wing analysis Is presented. 
Since the aerodynamic theory Is nonlinear, the 
method requires the coupling of two Iterative 
processes - an aerodynamic analysis and a 
structural analysis. A full potential analysis 
code, FL022, is combined with a linear 
structural analysis to yield aerodynamic load 
distributions on and deflections of elastic 
wings. This method was used to analyze an 
aeroelastlcally-scaled wind tunnel model of a 
proposed executive -jet transport wing and an 
aeroelastlc research wing. The results are 
compared with the corresponding rigid-wing 
analyses, and some effects of elasticity on the 
aerodynamic loading are noted. 

Nomenclature 

a local speed of sound 

b wing span 

c local chord 

C influence function; deflection at 

(x,y) due to a unit load at ( 5 , n) 
C zz deflection at y due to a unit load 

at n . 

C 60 angular rotation at y due to a unit 

moment at n 

c 1 section lift coefficient 

Cl wing lift coefficient 

Cm wing pitching moment coefficient 

Cp pressure coefficient 

aCp lifting pressure coefficient 

F lifting load element 

h wing deflection 

m twisting moment about elastic axis 

M Mach number 

q dynamic pressure 

u,v,w velocities, normalized with respect 

to free stream speed, in the x,y and 
z directions, respectively 
x,£ streamwlse coordinate 

x, 5 coordinate perpendicular to elastic 

axi s 

y, n spanwise coordinate 

coordinate along elastic axis 
z coordinate normal to x-y plane 

a wing root angle-of -attack 

Y ratio of specific heats 

e streamwise twist angle 

* velocity potential 

<d relaxation factor 

Introduction 

Methods for predicting transonic 
aerodynamic loads on wings have been developed 


assuming complete wing rigidity.! -3 in 
reality, airplane wings are flexible and, as Is 
commonly known, experience deflections that may 
significantly alter their loading. Thus, 
accurate prediction of aerodynamic loads on 
flexible wings and proper Interpretation of 
experimental data require consideration of the 
effects of elasticity on these loads. 

When the flow field Is described by linear 
equations, aeroelastlcally-corrected loads can 
be easily determined by direct solution of a set 
of matrix equations relating angle-of -attack, 
dynamic pressure, and structural and aerodynamic 
Influence coefficients. 4 For example. Pal and 
Sears® used matrix methods to calculate lift 
distributions on swept, flexible wings and 
demonstrated some effects of sweep on aero- 
elastlc phenomena. The FLEXSTAB program** 
represents the state-of-the-art in methods for 
computing aerodynamic loads on elastic wings. 

At transonic speeds, however, such direct 
methods are not yet developed. Alternative 
methods include correcting calculations made 
using matrix methods with empirical relation- 
ships derived from experimental transonic 
data? and using deflected wing shapes measured 
during wind tunnel tests in transonic analysis 
codes.® Since these methods require support- 
ing experimental data, they are quite limited In 
their application. 

Chipman et al.® developed a procedure for 
including elastic effects in transonic load pre- 
dictions that requires no experimental data. 

That procedure iterates between a transonic 
small disturbance method and a linear structural 
approximation that models wings as slender 
beams. Using not-fully-converged aerodynamic 
loads, deflections and rotations of streamwlse 
strips of wings are calculated and used to de- 
fine new wing shapes. Those shapes are then 
used in the aerodynamic computations, and the 
process is repeated until a converged aero- 
dynamic solution and wing shape are obtained. 

The choice of a beam structural model limits the 
method to the analysis of high-aspect-ratlo 
wings with chordwise rigidity. Also, as the 
authors note, for some supercritical Mach num- 
bers, small disturbance theory may predict shock 
waves at incorrect locations, causing moments 
and calculated twist angles to be in error. 

These points suggest the need for a more 
accurate aerodynamic theory and a more general 
structural representation. 

This paper presents an Improved method, 
based on full potential aerodynamics, for cal- 
culating steady, transonic loads on flexible 
wings. The present method also Iterates between 
a nonconverged aerodynamic calculation and a 
linear structural analysis. Since, the aero- 
dynamic analysis is based on full potential 
theory, calculated moments and related twists 
should be more accurate than those obtained 




1 



with small disturbance aerodynamic loads. 
Assuming linear structural relationships and 
hence small wing deflections, the structural 
analysis employs an Influence coefficient method 
which can be used without the slender beam 
assumptions. In computing deflections, 
distributed loads are represented as a network 
of discrete loads located at the centroids of 
corresponding area elements. Representing the 
loads In this manner allows wing deflections to 
be computed directly at the load points without 
the beam assumptions. The present method can, 
of course, be used for beam analysis as a 
particular case. To demonstrate the present 
method, aerodynamic loads on an 
aeroelastlcally-scaled wind tunnel model of a 
proposed executive-jet transport w1ngl° and on 
an aeroelastlc research wing* 1 have been 
computed. Since bending and torsional 
stiffnesses along the elastic axis were 
available for these wings, slender beam theory 
was used for the structural analyses of the 
wings studied in this paper. The results of 
those analyses are presented and compared with 
corresponding rigid-wing analyses and measured 
data. 


Aerodynamic Analysis 

Aerodynamic solutions are computed numeri- 
cally using the FL022 computer program 2 which 
solves the full potential equation 



-2(uv<f X y + uw* xz + vw$ yz ) ■= 0 (1) 

where a is the local speed of sound, 
and (u,v,w) = ($ x , <ty, $ 2 ) are the local 
components of velocity normalized with respect 
to the free stream speed. Equation (1) is 
solved, using a nonconservative finite 
difference procedure, on a sequence of 
successively finer grids. This procedure is 
computationally more efficient than using a 
single fine grid. 

For the results presented in this paper, 
the initial grid was chosen to have 48, 8, and 8 
points in the x, y, and z directions, 
respectively, and was eventually refined to 192 
x 32 x 32. The results presented herein were 
obtained on the finest grid. All aerodynamic 
computations were performed on a Control Data 
Corporation (CDC) CYBER 203 using a partially 
vectorized version of FL022. Nondimenslonal 
pressure distributions, C p , were then obtained 
from the following expression 


C = — {[ 1 + M 2 (l -(u 2 +v 2 +w 2 )) ] -1} 

yM (2) 

It should be noted that FL022 cannot be used to 
calculate the effects of fuselage interference 
on wing loading. 

Structural Analysis 

For a given load distribution, ap(x,y), 
wing deflections, h(x,y), can be obtained from 


h(x,y) ■ //C(x,y;e,n)APU.n)d5«ln (3) 


where the flexibility Influence function, 
C(x,y;S,n), defines the deflection at (x,y) due 
to a unit load at (?,n). Since the elastic pro- 
perties of a typical wing are generally very 
complicated. It Is usually not possible to 
derive analytic expressions for Influence func- 
tions. Thus, It Is usually necessary to deter- 
mine numerical values of C at specific locations 
on the wing. In the present analysis, load 
distributions obtained by solving (1) are 
represented as a network of load elements 
(discrete forces) concentrated at the centroids 
of the areas as shown In figure 1. Deflections, 
{h} , at those points can then be expressed as 

{h} - [ C ] { F} (4) 

where Fi = (ApA£dn)l is the force applied at 
the centroid of the 1 th area element. 

An Important advantage of representing the 
distributed load as a network of discrete forces 
is that it allows the structural analysis to be 
performed without the slender beam assumptions. 
That is, given an Influence coefficient matrix, 
deflections at the load points can be computed 
directly from (4). This means that low-aspect- 
ratio wings, which generally cannot be modeled 
as slender beams, may be analyzed using the pre- 
sent method. The structural influence coeffi- 
cient matrix, C, can be determined, for example, 
using a finite element structural analysis or by 
experiment. However, because linear structural 
relationships are used, the present analysis is 
limited to cases where the wing deflections are 
small. 

For the applications presented In this 
paper, the wings were of high aspect-ratio and 
In the present analysis were treated as slender 
beams. The beam analysis requires the defini- 
tion of an elastic axis and distributions along 
the elastic axis of stiffnesses In bending (El), 
torsion (GJ), and shearing (GK). Each load 
element is represented as a force at the elastic 
axis and a moment about the elastic axis (see 
figure 2). For slender beams, the influence 
function can be written as 


C(x ,y;C»n) = C zz (y,n) + xC 00 (y,n) (5) 

where C zz is the vertical deflection at y due 
to a unit load at rf, and C 09 is the angular 
rotation of a strip normal to the elastic axis 
at y due to a unit moment at n. C zz and 
C 00 are derived from the following expressions 
given in chapter two of reference 1 
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where X Is a dummy variable of Integration, and 
the coordinate system Is defined In figure 2. 

GK distributions were not available for the 
wings analyzed In this study, therefore shearing 
effects were not Included In C zz . Since (6) 
and (7) are used to derive C zz and C 00 along 
the elastic axis, the structural influence 
coefficient matrices can be chosen to have a 
desired number of elements which In the present 
analysis Is equal to the number of load ele- 
ments. Using the beam theory assumptions, 
deflections, { h} , at the locations of the load 
elements are given by 


{h} = [ C zz ] { F} + r<| [C 08 ]{m} (8) 


where mi = Si Fi Is the twisting moment 
about the elastic axis associated with each load 
element. Equation (8) Is valid only for wings 
that can be modeled as slender beams and thus is 
a special case of (4). 

The deflection calculation is performed 
external to FL022, on a CDC CYBER 175, and is 
Iterated as the aerodynamic load Is updated 
(typically, deflections are computed after each 
20-25 Iterations of the aerodynamic solution). 
Since the aerodynamic calculation is itself 
Iterative, the two Iterations are converged con- 
currently (see figure 3). It Is relatively 
straightforward to Include the deflection and 
aerodynamic load calculations In a single code, 
but the present procedure allows periodic exam- 
ination of the intermediate solution. The cost 
of a converged flexible-wing calculation Is not 
a great deal more than that of corresponding 
rigid-wing computations. This Is because a not- 
ful ly-converged aerodynamic solution is used for 
each deflection calculation and the computation 
of deflections requires relatively few computer 
resources. 

At the n th structural Iteration, inter- 
mediate deflections are computed from (8) and 
new deflections determined using the relaxation 
formul a® 

{h} n+1 = {h} n + co({h) -{h} n ) (9) 

where { h } represents the Intermediate deflec- 
tions. In this study the deflections were 
always underrelaxed (w=0.75). The relaxed 
deflections are added to the original wing 
coordinates to obtain a new wing shape and the 
Iterative solution of (1) by FL022 Is continued. 
As In reference 9, this process (figure 3) Is 
repeated until a converged aerodynamic solution 
and wing shape are obtained. In this study, 
solutions were considered converged when * 


changed by less than 2x10-5 between 
aerodynamic Iterations and all wing deflections 
changed by less than two percent between 
structural Iterations. 

The calculations were performed subject to 
one of the following conditions: (1) wing root 
angle-of -attack (a) specified and the resulting 
loads and deflections calculated or (2) wing 
lift coefficient (Cl) specified and the 
required root angle-of -attack determined during 
the calculations. In the latter case, because 
FL022 has no provision for specifying lift, 
necessary adjustments In angle-of-attack were 
made at each structural Interaction until the 
solution converged to the desired lift. 


Appl Icatlons 

The present method has been used to cal- 
culate loads and deflections for a wind tunnel 
flutter model of a proposed executive-jet- 
transport wing 10 and for the DAST (Drones for 
Aerodynamic and Structural Jesting) 1S[RW-2 (Aero- 
elastlc Researching Number J).U The 
planforms of these wings are shown In figure 4. 
For the transport wing model, some comparisons 
of calculated and measured tip deflections and 
twist angles are shown. While Influence 
coeffl cents were not measured for either wing, 

El and GJ distributions were available. 
Therefore, beam theory was used for the 
structural analyses of these high-aspect-ratio 
wings. 

For use In the deflection calcuatlons, 
pressure distributions obtained from FL022 were 
represented as a network of 200 discrete load 
elements with 10 elements per chord at 20 
locations on the semi span. Equation (8) was 
then used to compute deflections at each load 
element location. On the CYBER 175, two times 
as many load elements per span location could 
have been used, at little extra cost, with the 
calculations still accomplished In memory. 


Transport Wing Model 

The transport wing model, which had super- 
critical airfoil sections, was a 1/6.5 scale 
flutter model of a proposed executive-jet trans- 
port wing and was designed for Mach number (M) 
of 0.82 and dynamic pressure (q) of 30 psf. It 
was constructed of fiberglass front and rear 
spars and fiberglass ribs to which fiberglass 
skins were bonded. To prevent buckling, 
half-inch thick foam plastic panels were bonded 
to the interior of the skins between the ribs 
and spars. Measured bending and twist slopes 
about the elastic axis were used to determine El 
and GJ distributions. 


ARW-2 

The second DAST aeroelastlc research wing 
is of high-aspect-ratio (10.3) and has super- 
critical airfoil sections. Its planform Is 
similar to that considered optimum for an energy 
efficient transport. The structure consists of 
aluminum spars located at the 25 and 62 percent 
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chord lines, fiberglass skin panels riveted and 
bonded to the spars, and leading and trailing 
edges that are attached with screws. The wings 
are attached to the fuselage with a nearly rigid 
aluminum carry -through center section. A stress 
analysis has been performed to determine an 
elastic axis and El and GO distributions. 


Results 

Transport Wing Model 

Figure 5 shows the computed load distri- 
bution on the flexible transport wing model at 
five percent semispan Intervals for the design 
condition. The wing Is primarily aft-loaded, 
with small positive loading at the leading edge 
out to approximately the 70 percent semi span 
station. Outboard of this station, the load is 
negative near the leading edge and positive 
toward the trailing edge. Calculated deflec- 
tions along the elastic axis and streamwlse 
twist angles due to aeroelastic deflection, 
along with the tip deflections and twists 
observed during a wind tunnel test of the 
model, 12 are shown In figure 6. The calcu- 
lated twists were determined by dividing the 
differences in leading and trailing edge deflec- 
tion by the local chords. The experimental 
twists were measured optically with a catheto- 
meter. Agreement between the calculated deflec- 
tions and twist angles and those observed during 
the wind tunnel test Is very gbod. The observed 
tip twist Is only slightly larger than the 
calculated value. Contours of constant 
deflection, presented in figure 7, show that the 
wing deformation is composed primarily of 
bending of the elastic axis. It is well known 
that bending of a sweptback wing results in 
washout along the semispan i.e., negative 
Induced twist angles (figure 6b). Torsional 
deformations are relatively small and increase 
along the semi span as evidenced by the change in 
angle, from root to tip, between the deflection 
contours and the elastic axis. 

These results are very similar to those, 
for the same wing, presented in reference 9 
which were calculated with small disturbance 
aerodynamics. To place this comparison in 
proper context, however, it is observed that (1 ) 
small disturbance theory should yield its best 
results near design conditions where 
perturbations are not large and shocks are weak 
and (2) when using small disturbance theory, 
pressure differences, which cause structural 
deformations, are expected to be computed more 
accurately than actual surface pressures. At 
off -design conditions, where perturbations may 
be large, methods which use full potential 
aerodynamics should more accurately predict 
structural deformations and surface pressures. 

Figure 8 shows twist angles calculated with 
the present method at a high dynamic pressure. 
Also shown for comparison are the tip twist 
measured during the wind tunnel test, and tip 
twists presented In reference 9, which were 
obtained using small disturbance aerodynamics, 
with and without viscous corrections, and also 
with an Initial Inviscid full potential 
calculation. 9 in this case, the use of small 


disturbance aerodynamics, even with viscous 
corrections, results in calculated tip twists 
that are less accurate than those obtained with 
Inviscid full potential theory. The tip twist 
determined using the present analysis Is nearer 
to the measured value than that obtained by 
Chlpman et al. using full potential aerodynamic 
loads. Because of the limited size of the 
measured structural Influence coefficient matrix 
used In reference 9, wing deformations were com- 
puted at relatively few spanwlse locations and 
Interpolated to the desired locations. The 
closer agreement of the present result and 
experimental data Is possibly because the 
deflections were calculated directly at many 
more wing locations. Including viscous correc- 
tion In the small disturbance load calculations 
resulted In a more accurate predicted tip twist 
angle (figure 8). This suggests that at high q. 
Incorporating viscous effects In the present 
analysis should lead to a calculated tip twist 
that Is very near the measured' value. 

The calculated load distribution at the 
higher dynamic pressure Is shown In figure 9. 

At outboard stations, there are large negative 
loads on the forward part of the wing due to 
high suction pressures ahead of a strong shock 
near the lower surface leading edge. To Illus- 
trate this loading more clearly the surface 
pressure distribution at the 90 percent semispan 
station Is shown In figure 10. The high suction 
pressure on the lower surface Is quite promi- 
nent. 

The calculations required to analyze the 
transport wing are shown in Table I. Computa- 
tional requirements are listed In terms of the 
number of Iterations, on the fine grid, of the 
aerodynamic solution required to obtain a con- 
verged solution. While no rigid-wing results 
are presented In this paper, the number of iter- 
ations required to perform a rigid-wing analysis 
is listed to illustrate the additional computa- 
tions required for the flexible-wing analyses. 

The additional computations, 40 percent for the 
design case and 56 percent for the off -design 
case, are reasonable for the improved results 
obtained from the flexible wing analyses. 


ARW-2 

Transonic aerodynamic loads on the ARW-2 
wing were calculated at the design cruise condi- 
tions, M = 0.80, q = 126.4 psf. Comparisons of 
lifting pressures and actual surface pressures 
calculated assuming rigid and flexible wings at 
constant a and at the ARW-2 design cruise Ci 
are among the results that are presented. Since 
fuselage effects cannot be calculated using 
FL022, fuselage lift and aerodynamic Inter- 
ference were neglected in this analysis. The 
DAST ARW-2 configuration has a design cruise 
C|_ of 0.53. An estimated stabilizer trim 
angle of -4.2° and a stabilizer lift-curve 
slope of 0.012 were used to determine a cruise 
trim Cl of 0.58 for the wing. Because ARW-2 
is joined to the DAST vehicle with a nearly 
rigid carry-through center section, cambering of 
the wing due to fuselage deflection wa # s assumed 
to be negligible. Therefore, it was considered 



reasonable to analyze the ARW-2 wing In Its free 
flight trim condition as a flexible, canti- 
levered wing with Its root angle-of-attack 
chosen such that C|_ = 0.58. Wind tunnel data 
which was modified by a FLEXSTAB analysis was 
used to determine an estimated wing trim a of 
1.36°. In the present inviscld analysis, how- 
ever, the design trim C[_ was obtained at 
a = 0.91° for the flexible wing and a = 0° 
for the rigid wing. 

Figure 11 shows calculated rigid-wing and 
flexible-wing lift and pitching moment (Cm) 
coefficients for a range of angles-of-attack at 
M = 0.80, q = 126.4 psf. As expected, the 
inclusion of flexibility effects in the analysis 
results in considerable load relief on the 
wing. The calculated spanwise lift distri- 
butions shown In figure 12 Illustrate the loss 
of lift due to aeroelastic deformations at a = 
1.36°. In this case. Including aeroelastic 
effects In the wing analysis resulted in a 20 
percent reduction in the calculated lift. The 
effects of wing flexibility on the nondlmen- 
sional lifting pressure (AC p ) distributions Is 
shown In figure 13 at intervals of five percent 
semispan. The rigid-wing pressure distribution 
(figure 13a) Is characterized by a strong upper 
surface shock along the entire semi span, as 
evidenced by the humps in the ACn distri- 
butions, and relatively high lift near the lead- 
ing edge. When flexibility effects are Included 
in the analysis, the shock moves forward and 
weakens inboard of approximately the 50 percent 
semi span station and outboard of the 80 percent 
semi span station and completely disappears bet- 
ween the 50 percent and 80 percent stations 
(figure 13b). A comparison of rigid-wing and 
flexible-wing pressure distributions at the 90 
percent semi span location (figure 14) shows in 
detail the effects of aeroelastic deformations 
on surface pressures and on shock strength and 
location. Aft of 70 percent chord, however, the 
pressure distribution is virtually unchanged. 

Results of the analysis of ARW-2 at Cl = 
0.58 are presented In figures 15-18. Calculated 
spanwise lift distributions on the rigid wing at 
a = 0° and the flexible wing at a = 0.91° 

(figure 15) show how the lift is shifted inboard 
on the flexible wing. Figure 16 shows the load 
distributions at intervals of five percent semi - 
span. The more pronounced effects are seen at 
the outboard stations where, near the leading 
edge, the lift on the flexible wing is signifi- 
cantly less than that on the rigid wing. This 
is due to the smaller suction pressure at the 
upper surface leading edge caused by a decrease 
in the local angle-of-attack. A plot of surface 
pressures at the 90 percent station (figure 17) 
shows that aeroelastic effects cause the upper 
surface shock to weaken and move forward and the 
suction peak to decrease by about 18 percent. 
However, the pressure distribution on the aft 30 
percent of chord is unchanged. As in the case 
of the transport wing model, at the ARW-2 design 
Cl, bending of the elastic axis accounts for 
most of the wing deflection (figure 18). 

The number of fine-grid iterations of the 
aerodynamic solution required to analyze ARW-2 


are presented in Table II. The rigid-wing 
analysis at a « 1.36° required more fine grid 
iterations than the corresponding flexible wing 
analysis because very few calculations were made 
on the two coarser grids before doing the fine 
grid calculations. This slowed convergence of 
the fine grid solution. A more representative 
example Is the 100 fine grid iterations required 
for the rigid wing analysis at 0°, in which 
substantially more calculations were made on the 
coarser grids to Increase the rate of 
convergence. Using 100 as the number of fine 
grid Iterations required to obtain a rigid-wing 
solution, the flexible-wing solutions for 
a = 0.91° and 1.36° require 15 percent and 
35 percent more iterations, respectively. The 
additional cost of these calcuatlons Is quite 
reasonable. 


Concluding Remarks 

A method for calculating steady transonic 
loads on flexible wings has been developed by 
combining a nonlinear full potential flow 
analysis with a linear structural analysis. The 
structural model is not limited by the 
assumptions of beam theory. Because of the type 
of structural information that was available, 
however, the present method was used in a beam 
analysis mode to analyze two high-aspect-ratio 
swept wings - a wind tunnel flutter model of a 
proposed executive-jet transport wing and the 
DAST (Drones for Aerodynamic and Structural 
Jesting) ARW-2 (Aeroelastic Research Wing Number 
2). At design conditions where flow 
perturbations are relatively small and shock 
waves are weak, calculated twist angles and 
deflections of the tranport wing model show good 
agreement with measured data and with previously 
published results obtained using transonic small 
disturbance aerocfynami cs. At off-design 
conditions, however, where perturbations are 
larger and shocks are stronger, results 
determined using full potential aerodynamics 
were shown to be more accurate than those based 
on small disturbance theory. 

Results for the transport wing model at an 
off-design condition indicate that the ability 
to model the wing and its loading as a desired 
array of area elements and discrete loads, as in 
the present structural analysis, can lead to 
increased accuracy of the computed aeroelastic 
deformations. For the two wings studied in this 
investigation, aeroelastic effects on the load 
distributions at design conditions were 
primarily due to bending deflections. These 
bending deflections of swept wings induced 
negative streamwise twist angles and the 
resulting changes in calculated loads. 

The flexible wing analyses presented in 
this paper required from 15 to 56 percent more 
computations than the corresponding rigid-wing 
analyses. Even with the Increased computing 
requirements, the costs of the flexible wing 
analyses were still quite reasonable. 
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TABLE I. COMPUTATIONS REQUIRED TO ANALYZE 
TRANSPORT WING MODEL 




Aerodynamic 

Structure 

Flow Conditions 

Iterations 

Rigid 

Design Cruise 

90 

Flexible 

Design Cruise 

125 

FI exlble 

High q 

140 


TABLE II. COMPUTATIONS REQUIRED TO ANALYZE 
ARW-2 


Structure 

- 

Flow Conditions 

Aerodynamic 

Iterations 

Rigid 

Design Cruise, <*=1.36° 

190 

Flexible 

Design Cruise, a=1.36° 

135 

Rigid 

Design Cruise, a=0° 

100 

FI exlble 

Design Cruise, a=0.9l° 

115 
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^FREE STREAM 




Fig. 3 Aerodynamic/structural Iteration 
process. 



Fig. 2 Reference axes and representation of 
load element. 




(b) ARW-2 

Fig. 4 Planforms of wings analyzed in present 
study (dimensions in Inches). 
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(a) Lift coefficients 
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(a) Rigid wing 



STREAM 


( b ) Flexible wing 

Fig. 13 Nondimensional loading on ARW-2 at 
M - 0.80, a * 1.36°, q - 126.4 psf. 
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a, deg 

(b) Pitching moment coefficients 

Fig. 11 ARW-2 rigid-wing and flexible-wing lift 
and pitching moment coefficients at 
M = 0.80, q =126.4 psf. 
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Fig. 12 Spanwise lift distributions on rigid 
and flexible ARW-2 at M = 0.80, 
a = 1.36°, q = 126.4 psf. 


Fig. 14 Nondimensional pressures on rigid and 
flexible ARW-2 at ^ = 0.90 and at 
M = 0.80, a = 1.36°, q = 126.4 psf. 




O RIGID, fl = <r 
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\ /-RIGID 


_2l FLEX,BLE 


Fig. 15 Spanwlse lift distributions on rigid 
and flexible ARW-2 at M - 0.80, 

C[_ * 0.58, q = 126.4 psf. 
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(a) Rigid wing (a = 0°) 


STREAM 


0 .2 .4 .6 .8 1. 

x/c 

Fig. 17 Nondlmenslonal pressures on rigid and 
flexible ARW-2 at |£ = 0.90 and at 
M = 0.80, Cl = 0.58, q = 126.4 psf. 
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(b) Flexible wing (a = 0.91°) 

Fig. 16 Nondlmensional loading on ARW -2 at 
M = 0.80, Cl = 0.58, q = 126.4 psf. 


Fig. 18 Contours of constant deflection (In 
Inches) of ARW-2 at M = 0.82, 

Cl = 0.58, q = 126.4 psf. 
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